On the magnetic equation of state in (2+l)-flavor QCD 
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A first study of critical behavior in the vicinity of the chiral phase transition of (2+l)-flavor QCD 
is presented. We analyze the quark mass and volume dependence of the chiral condensate and chiral 
susceptibilities in QCD with two degenerate light quark masses and a strange quark. The strange 
Oh' quark mass (ms) is chosen close to its physical value; the two degenerate light quark masses (mi) are 

varied in a wide range 1/80 < nii/ms < 2/5, where the smallest light quark mass value corresponds 
to a pseudo-scalar Goldstone mass of about 75 MeV. All calculations are performed with staggered 
QQ ' fermions on lattices with temporal extent A''^ = 4. We show that numerical results are consistent 

I with 0{N) scaling in the chiral limit. We find that in the region of physical light quark mass 

values, rui/ma ~ 1/20, the temperature and quark mass dependence of the chiral condensate is 
already dominated by universal properties of QCD that are encoded in the scaling function for the 
' chiral order parameter, the magnetic equation of state. We also provide evidence for the influence 

of thermal fluctuations of Goldstone modes on the chiral condensate at finite temperature. At 
^ temperatures below, but close to the chiral phase transition at vanishing quark mass, this leads to a 

^ ' characteristic dependence of the light quark chiral condensate on the square root of the light quark 
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I. INTRODUCTION 



I Chiral symmetry and its spontaneous breaking in the vacuum are key ingredients to our understanding of the phase 
^\ • structure of strongly interacting matter at non-zero temperature and vanishing baryon chemical potential. In the 

' limit of massless quark flavors the QCD phase transition is controlled by the SUL{nf) x SUnirif) chiral symmetry. 
On , Quite general renormalization group arguments suggest [1] that QCD with three degenerate light quark flavors has a 

■ first order phase transition, whereas the 2-flavor theory is expected to have a second order phase transition. In the 
latter case the SUl{2) x SUii{2) chiral symmetry is isomorphic to 0(4) and the transition therefore is expected to 
belong to the same universality class as 3-dimensional, 0(4) symmetric spin models. Depending on the value of the 

. _ strange quark mass the QCD phase transition in the limit of vanishing light quark masses (up, down) may be first 

■ order or a continuous transition still belonging to the 3-dimcnsional, 0(4) universality class [1]. 

■ ■ ' While numerical calculations in 3-flavor QCD gave evidence for the existence of a first order transition, many of 
the details of the transition in 2- or (2+l)-flavor QCD with light up and down quarks arc still poorly constrained 
through lattice calculations. In particular, we do not know whether the chiral phase transition in (2-|-l)-flavor QCD 
is first or second order. An answer to this question is not only of academic interest; it also greatly influences our 
thinking about the phase diagram of QCD at non-zero baryon chemical potential [2]. The present analysis, although 
still performed on rather coarse lattices, is a first step towards answering this question. 

Attempts to verify the universal critical behavior associated with the QCD chiral phase transition in 2-flavor QCD 
have been made already before in calculations with staggered [3-9] and Wilson [10, 11] fermions. None of these lattice 
discretization schemes for the fermion sector of QCD preserve the full chiral 0(4) symmetry of the QCD Lagrangian. 
It therefore may not be too surprising that the early attempts to verify universal scaling properties of QCD were not 
too successful. In fact, at non-zero lattice spacing the Wilson fermion formulation does not preserve any continuous 
symmetry related to the chiral sector of QCD. The staggered formulation preserves at least an 0(2) symmetry at 
non-zero lattice spacing that gives rise to a single massless Goldstone mode in the chiral limit. Nevertheless, direct 
determinations of critical exponents within the staggered discretization scheme [3, 4, 7] did not deliver the expected 
0{N) results. Of course, due to the explicit breaking of 0(4) symmetry in the staggered formalism at non-vanishing 
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lattice spacing one would not have expected to be sensitive to 0(4) scaling. However, 0(4) and 0{2) critical exponents 
are quite similar and one thus might have hoped to observe at least some generic evidence for 0{N) scaling^. 

In the same spirit, the magnetic equation of state, i.e. the scaling of the chiral order parameter as function 
of reduced temperature and quark mass, has been analyzed subsequently in calculations with staggered [5, 6, 8] 
and Wilson fermions [10, 11]. The studies performed with Wilson fermions gave some indication for 0(4) scaling. 
These calculations, however, had been constrained to the high temperature, symmetry restored phase and had been 
performed with rather large values of the quark mass. They therefore did not allow to perform a test of scaling in 
the symmetry broken phase and could not contribute to the question of how Goldstone modes influence the scaling 
behavior at low temperatures. This contribution of Goldstone modes is a prominent feature of the magnetic equation 
of state of 0{N) symmetric theories which has been analyzed in detail in 0{N) symmetric spin models [12-16]. In the 
case of staggered fermions one might have hoped to find at least evidence for 0{2) scaling^. In fact, not only critical 
exponents, but also the 0{2) and 0(4) magnetic equations of state are quite similar. Deviations from the scaling 
function, however, turned out to be large in the low as well as high temperature regions and even in calculations 
on lattices with rather small lattice spacings [5, 6]. The missing evidence for 0{N) scaling also left room for an 
interpretation of the scaling behavior of 2-flavor QCD in terms of a first order phase transition [8]. 

We will present here results from calculations with staggered fermions. Contrary to most earlier studies of scaling 
properties these calculations have been performed with an action that suppresses cut-off efii'ects induced by a non- 
zero lattice spacing (a) in finite temperature calculations. Thermodynamic quantities are O(a^) improved. A first 
analysis of Goldstone effects with this action has been performed for rather large quark masses in Ref. [18]. Our 
calculations have been carried out with smaller than physical light quark masses so that the lightest Goldstone mode 
is a factor two lighter than the physical value of the pion mass. This will allow us to address the basic features of 
the thermodynamics induced by Goldstone modes and to gain some control over the universal features in the vicinity 
of the chiral phase transition temperature. We will present numerical evidence that at finite temperature, in the 
symmetry broken phase, the dominant quark mass dependence of the chiral condensate arises from fluctuations of the 
Goldstone modes that lead to a square root dependence of the condensate on the light quark masses. Such a behavior 
is expected in 3-dimensional theories with global 0{N) symmetry [19-22]. We will also show that universal scaling 
properties of the condensate are consistent with the 3-dimensional 0{N) universality class. Furthermore, an analysis 
of scaling violations, induced by the regular part of the QCD partition function, suggests that the crossover transition 
in QCD with physical quark masses is already strongly influenced by contributions arising from the singular universal 
part of the QCD partition function. At present we arc not sensitive to differences between 0(2) and 0(4) scaling. 
However, we point out that a combined analysis of scaling functions for the order parameter and its susceptibility 
should provide unambiguous results on the universality class of the chiral transition in QCD. 

This paper is organized as follows. In the next section we summarize universal properties of 3-dimensional 0(2) 
and 0(4) symmetric spin models and introduce notations. In Section HI we present our data on the quark mass and 
temperature dependence of chiral condensates in (2-|-l)-flavor QCD. The main results on the magnetic equation of 
state are discussed in Section IV. In Section V we give a brief account of properties of susceptibilities of the chiral 
order parameter. Section VI contains our conclusions. In Appendix A, for the readers' convenience we compile the 
asymptotic forms and the interpolations used for the 0(2) and 0(4) scaling functions, as adopted from [15, 16]. The 
numerical data which this paper is based on are summarized in Appendix B. 

II. 0{N) SYMMETRY BREAKING 

In the limit of vanishing light quark masses QCD is expected to undergo a phase transition at some critical 
temperature Tc at which chiral symmetry gets restored. The light quark chiral condensate, {ipil))i, will vanish at this 
temperature. Its quark mass and temperature dependence in the vicinity of the critical point, [T,mi) = (Tc,0), is 
controlled by a scaling function that arises from the singular part of the partition function. One way of analyzing the 
non-analytic structure of the QCD partition function, which has been pursued in the past, is to study the so-called 
magnetic equation of state. Before continuing our discussion of critical behavior in QCD we briefly summarize basic 
scaling relations using the conventional spin model notation, where the order parameter is denoted by M and the 



^ We are dealing here with numerical calculations of a cutoff theory whose Lagrangian has a global 0(2) symmetry. The relevant symmetry 
of QCD in the continuum limit, on the other hand, is expected to be 0(4). For many aspects of the discussion of critical behavior 
presented here, the distinction between 0{2) and 0(4) is of no importance. In these cases we will gcnerically talk about 0{N) symmetric 
models. 

^ Convincing evidence for 0(2) scaling has been found in calculations with a masslcss staggered fermion action to which an irrelevant 
chiral 4-fermion interaction has been added [17]. 
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FIG. 1: Scaling functions for the universality classes of three dimensional 0(2) and 0(4) models (left). Crosses show the 0(4) 
scaling function with an argument z = \.1z. The right hand figure shows the scaling function for the chiral susceptibility 
introduced in Eq. 7. 



symmetry breaking field is denoted by H. The critical behavior of 0(2) and 0(4) spin models in three dimensions 
has been analyzed extensively in the past. We will follow here closely the discussion given in [15]. 



A. Magnetic Equation of State 

In the vicinity of a critical point regular contributions to the partition functions become negligible and the universal 
critical behavior of the order parameter M of, e.g. 3-dimensional 0{N) spin models, is controlled by a scaling function 
fa that arises from the singular part of the logarithm of the partition function, 

M{t,h) = h^/'faiz) , (1) 

with z = t/h^/^^ and scaling variables t and h that are related to the temperature, T, and the symmetry breaking 
(magnetic) field, H, 

t = — ^ , /i = — . (2) 

<o Tc ho 

Here /3 and S are critical exponents characterizing the approach of the order parameter M to the critical point when 
one of the scaling variables is set to zero, 

M = {-tf , h = 0, t<0 (3) 
M = h^/^ , t = (4) 

These relations also fix the normalization of the scaling variables t and h, i.e. they define the normalization constants 
<o and Hq introduced in Eq. 2. Equivalently one can fix Iq and Hq through normalization conditions for the scaling 
function fa, 

/g(0) = 1 , lim 7^ = 1. (5) 

As it is of some relevance for our later discussion we note here that the normalization conditions for the scaling 
function, fc which fix the scale parameters Stud ho, refer to values of the scaling variable z, that are infinitely 
apart. 

For 0{N) symmetric spin models in three dimensions the scaling functions have been analyzed in much detail using 
Monte Carlo simulations and renormalization group techniques. In Fig. 1 we show results for the 0{2) [13] and 
0(4) [12, 14] scaling functions obtained by using the implicit parametrizations given in Ref. [15] and compiled in 
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N 




7 


S 


C2 


Zp 


2 


0.349 


1.319 


4.780 


0.592(10) 


1.56(10) 


4 


0.380 


1.453 


4.824 


0.666(6) 


1.33(5) 



TABLE I: Critical exponents /3, 7, 5 and the universal constant £2 for the three dimensional 0(2) universality class are taken 
from [15]; for 0(4) we use the data from [16]. The three critical exponents are related through 7 = P{S — 1). The last column 
gives the location of the maximum of the scaling function [15, 16]. 

Appendix A. There are clear differences between both scaling functions. We note, however, that the manifestation 
of the differences between the 0{2) and 0(4) scaling functions in the limited range jzj < 5 shown in Fig. 1 relies also 
on the normalization of these scaling functions in the limit z — > — 00. Without information on the scaling function at 
z —00, i.e. in a numerical study that only has access to a limited range of z-values, the scale parameters to and 
Hq could easily be adjusted to make the 0{2) and 0(4) scaling functions almost coincide. For [zj < 5 this is shown 
by the crosses in Fig. 1 which have been obtained by rescaling the argument of the 0(4) scaling function, z — > 1.2z, 
i.e. through a change of zq by 20% in this interval. This will lead to a violation of the normalization condition at 
z = —00 by a factor 1.2^ ~ 1.07, for the 0(4) scaling function. 

The above discussion makes it evident that a numerical analysis of the magnetic equation of state alone, in a scaling 
regime as large as jzj < 5, will not allow to distinguish 0(2) scaling from 0(4) scaling, unless the numerical accuracy 
is extraordinary high. Additional information will be needed to distinguish 0(2) from 0(4) scaling. This can be 
achieved through accurate control over the chiral limit, z = —00, or through the analysis of other scaling functions 
like the scaling function, /^(z), for the susceptibility of the order parameter, 

XMit,h) = 1^ = i-/,V5-i/^(z) , (6) 

/x(^) = J(/gW-^/gW) ■ (7) 

This scaling function is shown in Fig. 1 (right). It has a maximum at Zp which together with some critical exponents 
of the 0{N) models is given in Table I. It is obvious from this figure that a simple rescaling of the scale parameter z 
cannot transform an 0(2) scaling function into that for 0(4). 

While different 0{N) symmetric models are characterized by universal scaling functions, the scaling variables have 
to be normalized properly. The scale parameters to and ho are not universal. They depend on the 0(N) symmetric 
model under consideration, the definition of the scaling variables and also on the absolute normalization of the order 
parameter M. For instance, a rescaling of the order parameter by a constant factor, M — > bM, can be absorbed 
in a redefinition of the normalization constants to b^^/^to and ho — > b^^ho. This leaves the argument z of the 
scaling function and the scaling function itself unchanged, z — > z and M/h^/^ Mjh}-!^ . For a given definition of the 

symmetry breaking field H the combination zq = h}'^^^ jto therefore remains unchanged and is unique for any 0{N) 
symmetric model, i.e. its value is characteristic for that theory. It, for instance, characterizes the _ff-depcndence of 
the pseudo-critical line of transition temperatures, Tp{H), determined from the peak position of the order parameter 
susceptibility, 

Tp{H) - Tc ^ Zp j^i/ps j-g^ 

Tc Zq 

In Section IV we will determine the corresponding scaling relation for the pseudo-critical line of (2+l)-flavor QCD 
from an analysis of the magnetic equation of state. 

B. Contribution of Goldstone modes 

The spontaneous breaking of the continuous 0{N) symmetry at low temperature gives rise to massless Goldstone 
modes. The fluctuations of these light modes are reflected in the non-analytic dependence of the order parameter on 
the symmetry breaking variable, h [19]. In three dimensions this leads to [19, 21, 22] 



M(t, h) = M(i,0) + C2{t)Vh + 0{h) for alH < 



(9) 
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FIG. 2: The 0(2) scaling function, faiz) compared to the asymptotic form, /g'(z) and the leading order term, {—z)^. 

This leading correction to the temperature dependence of the order parameter, which arises from a non-vanishing 
explicit symmetry breaking (h > 0), is also reflected in the leading correction to the scaling function fa. For large, 
negative values of z one has, 

faiz) ~ f^iz) = +?2/3(-z)"W2) ^ for z ^ -oo . (10) 

As has been discussed also in [15], Eq. 10 can easily be obtained from the magnetic equation of state derived by 
Wallace and Zia [19]. The universal amplitude C2 is also given in Table I. In the following we will not make use 
of the scaling behavior of fa in the opposite limit, z — *■ +00. We note, however, that in this limit fa is controlled 
by the critical exponent 7 = f3(S — 1), fciz) ~ R^z~'^ . The universal parameter i?^ has been determined for three 
dimensional 0(2) [23] and 0(4) [24] universality classes. 

The asymptotic form of the 0{N) scaling function, f^{z), gives an excellent approximation to /g(^) in almost 
the entire low temperature regime, t < 0. This is evident from Fig. 2 where we compare fa to /q as well as to 
the leading order form, (--z)^. For z < —2 differences between fc and Jq are less than 2% for 0{2) and less than 
1% for 0(4). The leading /i-dependent correction that arises from the presence of Goldstone modes thus gives the 
dominant contribution to the order parameter in this regime. In order to establish universal critical behavior through 
an analysis of the 0{N) magnetic equation of state for QCD it therefore will be crucial to establish the influence of the 
Goldstone mode on the quark mass dependence of the chiral condensate and its derivative, the chiral susceptibility. 

III. CHIRAL SYMMETRY BREAKING IN (2+1) FLAVOR QCD 

A. Quark mass and volume dependence of the chiral condensate 

We discuss here our calculations performed for (2+l)-flavor QCD on lattices with size X Nr. We have fixed 
the temporal extent, Nt- = 4, and performed calculations for different spatial lattice sizes N„ = 8, 16 and 32 to 
control finite volume effects. All our calculations have been performed with a tree level improved gauge action and 
an improved staggered fermion action (p4-action), which eliminates 0(0^) discretization errors in thermodynamic 
observables at the tree level. The value of the bare strange quark mass in lattice units has been fixed to rhs = 0.065. 
In earlier calculations of the equation of state [25] and the transition temperature [26], performed with the same 
improved gauge and staggered fermion actions, it had been shown that in the present (small) range of gauge couplings 
this value of the strange bare quark mass yields almost physical values for the masses of the strange pseudo-scalar 
meson and the kaon. The light quark masses have been varied in the range 0.0008125 < ifii < 0.026. The smallest 
value, mi/ms = 1/80, corresponds to a pseudo-scalar Goldstone mass of about 75 MeV. For each of the 6 quark 
mass values chosen in the above interval we have performed calculations at several values of the gauge coupling in 
the range 3.28 < (3 < 3.33. As will become clear later this covers a temperature range 0.96<T/Tc<1.06, with Tc 
denoting the transition temperature in the chiral limit on a lattice with fixed temporal extent. All calculations have 
been performed using the Rational Hybrid Monte Carlo (RHMC) algorithm. In most cases we collected 15.000 to 
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FIG. 3: The light quark chiral condensate in lattice units versus the ratio of the light and strange quark masses (left) and its 
square root (right). The condensates have been calculated on lattice of size x 4 at the values of the gauge couplings shown 
in the right hand figure. The spatial lattice size for the two lightest quark mass values, mi/rris = 1/80 and 1/40, is Ncr = 32, for 
mi/nis — 1/20 it is = 32 for (3 = 3.28. In all other cases {mi/nis — 1/10, 1/5 and 2/5) the spatial lattice size is Na — 16. 



40.000 trajectories with length of half a time unit. We give more details on our simulation parameters, the statistics 
collected and expectation values of chiral condensates in Appendix B. 

Whenever we convert results to physical units we use the scale setting and meson mass calculations performed in 
connection with our calculation of the equation of state [25] and the transition temperature [26] . The most central 
anchor point for our current analysis is the determination of the pseudo-scalar mass (wps) and the Sommer scale 
parameter rg for light quark masses mi/nis = 1/20 at a gauge coupling (3 = 3.30. This value of the gauge coupling is 
close to the critical temperature in the chiral limit and the light to strange quark mass ratio is close to the physical 
mass ratio. For this parameter set we find mp^a = 0.1888(6) and rg/a — 1.8915(59) in lattice units. This converts to 
rups = 150.2(3) MeV when using tq = 0.469 fm [27] as it has been done also in our earlier calculations. 

The main part of our analysis is based on calculations of the light and strange quark chiral condensates, 

1 1 d\nZ 1 1 , 
,T,^ 1 1 a In Z 1 1 , 1 , 

- 4iv|Arw = i^^^^^">' ^''^ 

where Di and Ds denote the fermion matrices for light and strange quarks, respectively. 

A first overview on our data sample is given in the left hand part of Fig. 3. This figure shows results for the light 
quark chiral condensate in lattice units, calculated for different values of the gauge coupling and 6 different values 
of the light quark mass. For each parameter set we only show results from the largest spatial lattices available. We 
comment on finite volume effects below. We also note that we will conclude later that the chiral phase transition 
temperature at m; = corresponds to f3c — 3.30. As is evident from Fig. 3 the light quark chiral condensate shows 
a strong quark mass dependence that is not consistent with a linear dependence on mi. In fact, at low temperatures 
the dominant quark mass correction seems to be proportional to \fmi. This is highlighted in the right hand part of 
Fig. 3. Results for the light and strange quark condensates are summarized in Appendix B. 

In order to make sure that the drop in (V'V')i, seen for small values of the quark mass, is not due to a too small 
spatial volume, we analyzed the volume dependence of our results by performing calculations on lattices with spatial 
extent N^, = 8, 16 and 32. We show some results from this analysis in Fig. 4. As expected, the volume dependence of 
the chiral condensate increases with decreasing value of the quark mass. We have, however, no evidence for a strong 
increase of the volume dependence close to the phase transition temperature (/3 ~ 3.3). This suggests that also for 
the smallest quark masses used, our results obtained on lattices of spatial size 32'^ are close to the infinite volume 
limit. We also note that the smallest quark mass value used on our largest lattices corresponds to mpsN„ ~ 3, where 
nips denotes the lightest pseudo-scalar meson mass, i.e. the Goldstonc meson in the staggered fermion formulation of 
(2-hl)-flavor QCD. 



7 



<W>i 


m|/ms=1/40 


(3=3.28 ^ 

3.29 (shifted +0.1) — • 

3.30 (shifted +0.3) ^ 


^ 1/VT^ " 



i.OO 0.02 0.04 0.06 0.08 0.10 0.12 0.14 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 0. 



FIG. 4: The light quark chiral condensate in lattice units versus the inverse volume, V = A^^i for two values of the gauge coupling 
in the low temperature phase, /3 = 3.28, 3.29, and close to the transition temperature, /3 — 3.30, for mi/nis = l/40(left). The 
right hand figure shows results for the three smallest values of the light quark mass at /? = 3.30, i.e. close to the chiral transition 
temperature. For these quark mass values the largest lattice size, Na = 32, corresponds to 3 < nipsNcr < 6. In the left hand 
figure data for {4>tp)i have been shifted by a constant as indicated in the figure. 



B. The chiral order parameter 



The chiral condensate, as introduced in Eq. 11, is an order parameter for the chiral phase transition. At non- 
vanishing light quark mass an additive and multiplicative renormalization is needed to define an order parameter 
in the continuum limit. In Rcf. [25] we introduced an order parameter where quadratically divergent, additive 
contributions, which are proportional to the quark mass, have been removed by subtracting a suitable fraction of 
the strange quark chiral condensate from the light quark condensate. We will use this observable also here^. To 
take into account the anomalous dimensions of the chiral condensate we introduce a multiplicative renormalization, 
using the strange quark mass. A similar procedure has been suggested for 2-flavor QCD [5]. Furthermore, we express 
this product in units of to make it dimensionlcss. We thus introduce as an order parameter for chiral symmetry 
restoration in (2+l)-flavor QCD 

M = nis - —{4>^)s) ■ (12) 

V Tils J 

For our scaling analysis at fixed Nt a renormalization of the order parameter is not at all necessary. We may as well 
analyze the scaling behavior of the non-subtracted chiral condensate introduced in Eq. 11. To check the consistency 
of our analysis, we will do so and in the following also utilize the non-subtracted order parameter 

Mb = N^ms{^P^)i , (13) 

where we have introduced the constant multiplicative factor, N^rhs, such that Mb and M agree in the chiral limit. 



IV. THE MAGNETIC EQUATION OF STATE 



Having introduced our numerical results for the light and strange quark chiral condensates and the subtracted 
(M) and non-subtracted (Mb) order parameters, we are now ready to discuss critical behavior in the vicinity of the 
transition temperature in terms of the magnetic equation of state. 



^ This does not remove divergencies that are logarithmic in the cut-off. In the free field limit these divergencies are proportional to 
{mi/T)^ and are therefore expected to be numerically small 



8 





(7Tl;/ms)max 


ho 


to 


Tc [MeV] 


zo 


M 


1/20 


0.0048(5) 


0.0048(2) 


195.6(4) 


8.5 (7) 


Mb 


1/20 


0.0022(3) 


0.0037(2) 


194.5(4) 


6.8 (5) 


M 


1/40 


0.0042(6) 


0.0047(2) 


195.3(4) 


8.0 (8) 


Mt 


1/40 


0.0025(5) 


0.0040(2) 


194.8(4) 


7.0 (6) 



TABLE II: Fit results for the scale parameters ho and to and the chiral transition temperature Tc using the 0(2) scaling 
function. The last column shows the combination of scale parameters zo — h^^^/to- 

A. Scaling analysis 

In the vicinity of the chiral phase transition temperature corresponding to vanishing light quark masses and for 
sufficiently small explicit symmetry breaking, the order parameter is expected to scale according to Eq. 1. We introduce 
the reduced temperature and external field variables t and h as in Eq. 2. The definition of t obviously carries over from 
the spin model context to QCD. For the relation between the lattice gauge coupling (3 and the temperature we exploit 
the parametrization of the Sommcr scale parameter tq that has been determined in our calculations for the equation 
of state [25] at mi/nis = 1/10. This takes into account deviations from asymptotic scaling of the QCD /3-function for 
the range of couplings used here. We checked that the entire scaling analysis presented here only for a small range 
of the lattice cut-off is not really sensitive to these corrections and could as well have been performed by determining 
t = {T — Tc)/Tc using the asymptotic 2-loop /3-function or other approaches followed in earlier studies [3, 6, 8]. 
Likewise this analysis is, of course, independent of any physical value used to set the absolute scale for rg. 

The symmetry breaking external field H is proportional to the light quark mass. Also here we take care of the 
anomalous scaling dimension of quark masses and express the light quark mass in units of the strange quark mass, i.e. 
we introduce H = mi/m^. An alternative, yet similar way to deal with the anomalous dimensions has been suggested 
in [5]. 

In the chiral limit, at finite value of the cut-off (fixed Nr) we expect the phase transition in (2-|-l)-flavor QCD to 
be either first order or to belong to the universality class of three dimensional 0{2) models. In our current analysis 
we did not find any indication for a strong volume dependence or meta-stabilities in the time evolution of the chiral 
condensates or other obscrvables. In particular, as is also evident from Fig 10 shown in Section V, we have no evidence 
for a strong increase in the chiral susceptibility as function of volume. Although we can, at present, not rule out a 
weak first order phase transition at smaller quark masses, we have no indications for that to happen. We therefore 
will compare our data on the order parameter to the universal 0(2) scaling function, i.e. the magnetic equation of 
state introduced in Eq. 1 with critical exponents [3 and S given in Table I. We start by determining the three free 
parameters to, ho and the transition temperature, from fits to the order parameter M. For this we use the three 
lightest quark mass values, rrii/nis < 1/20, leaving out the lowest temperature value, corresponding to /3 = 3.28, and 
the two highest temperature values, corresponding to /? = 3.32 and 3.33. For mi/iris = 1/80 and 1/40 these data are 
from lattices of size 32"^ x 4 while the mi/nis = 1/20 data set is taken from calculation on 16'^ x 4 lattices. 

A posteriori we find that this temperature interval corresponds to 0.97 < T/Tc < 1.03. In this small temperature 
interval and for the small quark mass regime, which docs include the light quark mass value corresponding to the 
physical pion mass, we find good agreement between the rescalcd order parameter and the 0{2) scaling function. 
This is shown in the left hand part of Fig. 5. The fit yields /3c = 3.300(1) for the critical coupling corresponding to a 
phase transition temperature'' Tc = 195.6(4) MeV. We obtain an equally good agreement with the 0(2) scaling curve 
from an analysis of the non-subtracted order parameter M^. This is shown in the right hand part of Fig. 5. A similar 
analysis using the 0(4) scaling function and critical exponents yields scale parameters that are similar to those of the 
0(2) analysis, the largest differences occurring for to which comes out to be about 20% larger. 

The scaling function /g(z) is defined in the limit t ^ 0, h ^ 0, keeping z = t/h^/^^ fixed. In this limit the two 
order parameters M and coincide. From our scaling analysis we therefore should find identical results for the scale 
parameters, i.e. the critical temperature Tc as well as the normalization constants to, ho, if this analysis has been 
performed sufficiently close to the chiral limit. We have performed the scaling analysis for two different cuts on the 



This value is in excellent agreement with our earlier analysis performed with the p4 action on lattices with temporal extent Nt = 4. We 
stress, however, that this transition temperature is not extrapolated to the continuum limit and, in fact, within the staggered fermion 
approach, the chiral extrapolation should be performed after the continuum extrapolation to recover eventually the anticipated 0(4) 
scaling behavior. 
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FIG. 5: Fit of the 0(2) scaling function to numerical results for the subtracted order parameter M (left) and the non-subtracted 
light quark condensate Mb (right). This analysis has been performed for results obtained in calculations with light quark masses 
mi/ms < 1/20 and gauge couplings in the interval (3 £ [3.285,3.31]. 



ratio of the light to strange quark masses, mi/nis < 1/20 and mi/iris < 1/40. The fit parameters obtained in these 
two cases from an analysis of data for M and Mb are summarized in Table II. We note that results for Tc and to are 
within errors independent on the cut on nii/ms and the observable used. The scale parameter Hq is more sensitive on 
the choice of order parameter. However, there is a tendency that results for ho obtained from M and Mb converge to 
a common value if the cut on mi/rrig is reduced. 

The symmetry breaking field introduced in Eq. 2 is given in terms of the ratio of light to strange quark masses. To 
compare our result for scaling functions of QCD with other (model) calculations it may be more convenient to express 
H in terms of meson masses. In the present quark mass and gauge coupling range we find the approximate relation 
H = mi/nis ~ 0.52 {mps/rriK)'^ ■ We therefore may write the scaling variable z as 

^ — (^)'fe)"'- <»' 

As discussed in Section II this allows to determine the scaling behavior of the pseudo-critical line determined by the 
peak in the scaling function of the chiral susceptibility, f^^, 

Tpinips) - ^ p / rrv y^'"" ^^^^ 

Tc zo \mK J 

Using Zp from Table I and the values for zq given in Table II we find 0.682p/zo ~ 0.1 — 0.2. These values, which are 
consistent with earlier determinations of the slope of the pseudo-critical line [26, 28], emphasizes the weak dependence 
of the pseudo-critical temperature on the pseudo-scalar meson mass. 

We stress that this analysis has been performed in QCD at one non-vanishing lattice spacing, i.e. in the cut-off 
theory. The cut-off dependence of the normalization constants, the scale invariant ratio zq and the subtle continuum 
limit need to be studied in the future. We emphasize again, however, that the above combination of normalization 
constants for the scaling variables is an invariant of QCD and depends, in the continuum limit, only on the strange 
quark mass value. 



B. Scaling violations 



The scaling behavior observed for the chiral order parameters analyzed in the previous section is, of course, expected 
to hold exactly only in the limit t ^ and /i ^ 0, keeping the ratio z = tjh}!^^ fixed. At non-zero values of t and h we 
expect to observe scaling violations that may arise from sub-leading corrections to the scaling function as well as from 
the regular part of the QCD partition function. These corrections also depend on the definition of the order parameter. 
In particular, the two order parameters, M and A/t, introduced here differ in the treatment of contributions that are 
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FIG. 6: The order parameters M (left) and Mj, (right) for all quark mass values, mijvris < 0.4, and all values of the gauge 
coupling, /3 £ [3.28, 3.33], used in this study. The scaling variables t and h used to compare with the 0(2) scaling function are 
taken from the fit to the light quark mass results shown in Fig. 5. 



linear in the light quark mass. In our analysis of the order parameter, performed in a larger temperature and quark 
mass interval, we clearly see these differences and their role in contributing to violations of scaling. This is shown in 
Fig. 6. Most prominent are effects arising from a too large quark mass value. These effects show up in the scaling 
plot as deviations from the scaling function in the region of small z, i.e. for large quark masses at fixed t. They 
lead to the sizeable displacement of results obtained for too heavy quarks from the scaling curve. Effects that arise 
because the temperatures chosen are too far away from the critical point, t = 0, are typically not that drastic in our 
data sample. We fitted the scaling violations to an ansatz 

A/(t, h) = h}l^ jG{tlh}l^^) + atth + 6i/i + 63/1^ + 65/1^ . (16) 

We also considered including a term quadratic in the reduced temperature t^K). This correction, however, turned 
out to vanish within the errors of our fits. 

The fits of both order parameters performed with the ansatz given in Eq. 16 are shown in Fig. 7. As expected, we 
find that corrections linear in to;/to<, are eliminated in M. The corresponding fit parameter h\ is zero within errors 
and we therefore have fixed it to be zero in the fit shown in Fig. 7 (left). For the non-subtracted order parameter A/f, 
this term gives the dominant finite quark mass corrections. Here we find h\ = 0.0013(3). 



C. Scaling of the chiral condensate 



We have seen in the previous section that order parameters constructed from the chiral condensate are well described 
by the magnetic equation of state for small enough values of the light quark masses, 'mi/m,s<l/20. We want to 
underscore this point here by displaying the order parameters not in their scaling form, but as a function of temperature 
in units of the transition temperature determined in the previous section. This is shown in Fig. 8. The curves drawn 
in this figure are taken from the scaling fits to the subtracted and non-subtracted order parameters shown in Fig. 5. 
They had been obtained from the numerical results for M (left) and Aff, (right) in the range mi /iris < 1/20 and 
T/Tc = 1 ± 0.03. 



D. Comparison with earlier calculations in 2-flavor QCD 



As mentioned in the Introduction, there have been earlier attempts to compare the quark mass and temperature 
dependence of the chiral order parameter with 0{N) scaling functions on lattices with temporal extent Nr — 4 
[6, 8, 9]. These calculations had been performed for 2-flavor QCD using unimproved gauge and staggered fermion 
actions. In Ref. [6] calculations with three quark mass values had been performed, m = 0.008, 0.0125 and 0.025. The 
last two masses are similar to the two mass values used in Ref. [8], i.e. m = 0.01335 and 0.0267. In fact, results for 
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FIG. 7: The 0(2) magnetic equation of state compared to results for the subtracted order parameter M (left) and the non- 
subtracted chiral condensate, Mb for light quark masses mi/nis > 1/10. Curves show fits to data at fixed mi/ms using the 
ansatz for scaling violations given in Eq. 16. Same symbols correspond to same values of the gauge coupling. 
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FIG. 8: The subtracted chiral order parameter, M, defined in Eq. 12, compared to the fit result for the magnetic equation of 
state (left). The right hand figure shows results for the un-subtracted, but normalized chiral condensate Mb defined in Eq. 13. 

chiral condensates obtained in these two calculations are in good agreement with each other. This also is true for 
calculations performed in [4] where rh = 0.02 has been used. All these calculations have been performed at values of 
the gauge coupling in the vicinity of the cross-over at the corresponding quark mass values. They therefore mostly 
explored the region of z > 0. 

In Fig. 9 we compare results for the chiral condensate obtained in 2-flavor calculations with unimproved gauge and 
fermion actions with our results obtained in (2-|-l)-flavor QCD with 0{a^) improved gauge and fermion actions. In 
this figure we use a log- log plot as has been done also in Ref . [6] . In the 2-flavor case the symmetry breaking field has 
usually been chosen as H = rhN-r while for the reduced temperature variable we used (T — Tc)/Tc = R{(3c) / RiP) — 1, 
with f3c — 5.2435 as estimate for the critical value of the gauge coupling in the chiral limit [8] and R(J3) denoting 
the 2-loop /^-function for 2-flavor QCD. In the log-log plot differences in the scale parameters ho and zq correspond 
to shifts in vertical and horizontal directions, respectively. We made no effort to optimize the choice of these scale 
parameters for the 2-flavor data set. In Fig. 9 we have positioned the data such that the crossover region roughly 
corresponds to the location of the maximum in the 0(2) scaling function /^(zp), with Zp ~ 1.56 (see also [9]); this 
required the choice zq ~ 12. 

When comparing results obtained with standard and improved gauge actions the difference in the shape of the data 
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FIG. 9: Scaling plot for the chiral condensate calculated with an improved staggered action in (2+l)-flavor QCD (this work) 
and the standard staggered action in 2-flavor QCD [6, 8]. The results are shown in a log-log plot. For the (2-|-l)-flavor data 
set labels indicate the ratio mi/nia, in the 2-flavor case we give the bare quark masses m. Data for the lightest quark mass are 
from [6]. Data for the other two quark mass values are from [8]. For further discussion see text. 

sets clearly is the most striking feature. Apparently the (2-|-l)-flavor data set is in good agreement with 0{N) scaling 
while the results obtained with the standard staggered action deviate strongly. Furthermore, there is no tendency for 
better agreement with decreasing quark mass. Closer to the continuum limit, i.e. for larger N^-, results obtained with 
the unimproved actions have a similar shape [6] but seem to get somewhat closer to the 0{N) scaling curve. 

To compare the quark mass values used in the 2-flavor QCD calculations with those of the present (2-|-l)-flavor 
study we note that for standard staggered fermions at gauge couplings close to the cross-over a bare quark mass 
m = 0.025 corresponds to a pseudo-scalar Goldstone mass irips ~ 350 MeV [30]. The lightest quark mass used in 
these calculations, m = 0.008, therefore corresponds to nips ~ 200 MeV, which is similar to the pseudo-scalar mass 
obtained in (2-|-l)-flavor QCD calculations for the light to strange quark mass ratio mi/nis = 1/10. At this value for 
the light quark mass we observe only mild violations of scaling in calculations with the improved gauge and fermion 
actions. 

The number of flavors as well as the quark masses arc different in the data sets compared in Fig. 9. Nevertheless, 
it seems unlikely that this is the origin of the observed differences. It appears more probable that cut-off effects in 
calculations with unimproved gauge and fermion actions cause the differences. 

V. SUSCEPTIBILITIES 

In Eq. 6 we introduced the susceptibility xm as the derivative of the order parameter M with respect to the 
symmetry breaking field. Its temperature and quark mass dependence is controlled by the scaling function fx{z) 
defined in Eq. 7, which is shown in the right hand part of Fig. 1 for 0(2) and 0(4). Similarly we can, of course, 
introduce the susceptibility XMb as a derivative of the order parameter Mb with respect to H . Taking derivatives 
with respect io h = mi/{mshQ), rather than with respect to the ratio of quark masses, mi/ms, obviously requires 
knowledge of the scale parameter hg which we have determined in the previous section. 

In the analysis of QCD thermodynamics on the lattice it is more customary to calculate light (xD and strange (Xm) 
quark chiral susceptibilities, which are defined as derivatives of the corresponding chiral condensates with respect to 
mi/T and nis/T, respectively 



(17) 
(18) 
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To construct the susceptibility xm "^c will also need to take into account a mixed chiral susceptibility, 



The susceptibility of the subtracted order parameter, M , is then obtained as 

dM 

= XMb - hoN^ms ms - hoN^ml ^ ^ , (20) 

where in the last equality we introduced the susceptibility XMb, of the non-subtracted order parameter Mb. With this 
we can construct the scaling function for the chiral susceptibility, 

/x(^) = XMh,/h'/'-' ^ hi" f ^) ' XM , (21) 

from xm and similarly also from XMh- 

The scaling functions constructed from either xm or XMb, differ by terms that vanish in the chiral limit at T^. 
These terms therefore characterize once more systematic differences that arise in the construction of scaling functions 
due to the presence of regular terms that vanish, once the appropriate scaling limits are taken. We show scaling 
hmctions constructed from both order parameters in Fig. 10. We stress that all parameters (Tci a-iid /iq) have been 
determined in our analysis of the order parameters themselves. No fits are therefore involved in the comparison of 
the 0(N) scaling functions with the numerical results for susceptibilities shown in this figure. 

It is obvious from Fig. 10 that violations of scaling are significantly larger for susceptibilities than for the order 
parameters. Susceptibilities extracted from M and Mf, still differ for mi/nis = 1/20 but start to become compatible 
within errors for mi/ms < 1/40. For mi/ms = 1/40 we also show results from calculations on two different lattice 
size, 32'^ X 4 and 16^ x 4. If any, the volume dependence of susceptibilities is small at this value of the quark mass. 

There is yet another difference between order parameter susceptibilities derived in QCD, where the order parameter 
is a composed operator constructed from fermionic fields, and 0{N) spin models, where the order parameter is the 
expectation value of a scalar boson field. The order parameter susceptibilities in QCD receive two contributions, 
usually called the disconnected and connected part of the susceptibility, 



with 

Xdis 
Xcon 



xl^2xT^ + xr, (22) 

^{((Tri.-)V(Tri.rV}, (23) 
-J2{Di\x,Q)D;\Q,x)) . (24) 



While the first term, the disconnected part of the light quark susceptibility, describes fluctuations of the light quark 
condensate and has a direct analogy in the fluctuations of the order parameter in an 0{N) spin model, the second term 
(Xcon) arises from the explicit quark mass dependence of the order parameter, the chiral condensate. The connected 
part is an integral over the (quark-line connected) correlation function of the (iso-vector) scalar operator, ip^. The 
integral has a rather subtle quark mass dependence. Since 6 > 2, however, Xcon/h^"^^ will vanish in the chiral limit. 
In this limit the connected part of the susceptibility will therefore not contribute to the scaling function which in 
turn will entirely be determined through the disconnected part of the light quark susceptibility. 

At non-vanishing values of the light quark mass, however, the non-vanishing connected part of the chiral suscepti- 
bility is responsible for additional scaling violations in In fact, the scaling violations due to the connected part 
are distinctively different between 0{2) and 0(4) symmetric theories [22]. It is only in the latter that fluctuations 
of Goldstone modes do not contribute to Xcon- In the case of 0{2) symmetric models Xcon is expected to diverge 
proportional to \/Vh just like the total order parameter susceptibilities will do. 

In the staggered formulation of QCD with 2 light quark flavors the lack of 0(4) symmetry in the Lagrangian is due 
to explicit symmetry breaking terms (taste violations) that disappear only in the continuum limit. Corresponding 
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FIG. 10: Susceptibilities constructed from the subtracted order parameter M (left) and the non-subtracted light quark chiral 
condensate Mb (right). The data give results from calculations in (2-|-l)-flavor QCD on lattices with temporal extent A'^^ = 4 
and light quark mass values mi/nis < 1/20 in the interval /3 G [3.28, 3.33]. For the mi/ms = 1/40 data sample we show results 
for two different spatial lattice sizes. Filled symbols correspond to Na — 32 and open symbols are for = 16. 



to the 0(2) spin models, at finite values of the cut-off the divergence of Xcon ~ I/aA^ in the chiral limit can thus 
be understood in terms of taste violating contributions to the scalar correlation function [31]. We will discuss these 
subtle aspects of susceptibilities of the order parameter, the influence of taste violating terms in the staggered action 
on scaling properties of these susceptibilities and the resulting cut-off dependence of in more detail in a forthcoming 
publication [32]. 



VI. CONCLUSIONS 



We have performed a new analysis of scaling properties of the light quark chiral condensate in (2 + l)-flavor QCD. 
We found that at fixed non-zero lattice spacing the chiral condensate calculated with improved staggered fcrmions 
shows scaling behavior in the chiral limit that is consistent with 0(2) scaling. 

Through the analysis of scaling properties with quark masses that are smaller than the physical light quark masses 
we could fix the normalization constants to and Hq in the scaling variables t and h. This allowed us to quantify scaling 
violations for non-zero values of the quark masses in the vicinity of the phase transition temperature. These scaling 
violations turned out to be small in the magnetic equation of state already for physical values of the quark mass. 

On the basis of studying just the magnetic equation of state, we gave arguments that it will remain difficult to 
rule out 0(4) scaling without extraordinary precision of numerical lattice data. However, a distinction between 0(2) 
and 0(4) scaling might become possible through an accurate analysis of susceptibilities of the order parameter. At 
present, we still find significant deviations from scaling for the chiral susceptibility. This will be discussed in more 
detail in a forthcoming publication [32]. 

A determination of to and ho also fixes the scale parameter zq = hy^^ /to, which controls the quark mass dependence 
of the pseudo-critical line determined from the peak in the chiral susceptibility. In our present analysis this parameter, 
which uniquely characterizes non-universal aspects of critical behavior in QCD. has only been determined at one value 
of the lattice cut-off. Calculations at smaller lattice spacings. together with good control over scaling violations induced 
at non- vanishing quark masses will be needed to extract zq in the 0(4) symmetric continuum limit. 

The good scaling properties found here in calculations with 0{o?) improved gauge and fcrmion actions are in 
contrast to earlier calculations that had been performed with unimproved staggered fcrmion and gauge actions. We 
argued that the observed differences are due to cut-off effects. 

In our analysis we have assumed that the strange quark mass in (2 + l)-flavor QCD is large enough to avoid a 
first order phase transition in the light quark chiral limit. Although the good scaling properties of the chiral order 
parameters and the absence of a strong volume dependence in the light quark susceptibilities support this assumption, 
we clearly cannot exclude a first order transition to occur at still lighter quark masses. Consistent with limits given 
on the location of a first order transition in 3-flavor QCD [33, 34], however, our current analysis rules out such a 
transition for pseudo-scalar masses rups > 75 MeV. 
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APPENDIX A: SCALING FUNCTIONS FOR THREE-DIMENSIONAL 0(2) AND 0(4) MODELS 



In this appendix we summarize the scaling functions for models in the three-dimensional 0(2) and 0(4) universality 
classes. These interpolating functions have been taken from Ref. [15] and [16]. We note that the original parameters 
for the 0(4) model published in the tables of Ref. [15] had been updated in [16]. Moreover, interpolating curves had 
been constructed only for the scaling function fa- Applying these interpolations also to required slight adjustments 
of the interpolation parameters (j/Oj p)- 

In Eq. 1 we expressed the dependence of the order parameter M on the scaling variables t and h in terms of a 
scaling function, fa- Following the discussion given in [15] we introduce the variables x and y, 

y^fc': x^{t/h'/^')f^'/r (Al) 

Obviously y > 0. For small and large values of y the asymptotic forms that relate a; to y are known. For small y we 
have 

xM = -^ + (ci +d3)y + C2 + y'^' , (A2) 

and for large values of y one finds 

x,(?/) =ayi/^ + 6y(i-2/3)/7 , (A3) 
One can smoothly interpolate between these two relations [15] using the ansatz 



x{y) = Xsiy) 



P 



— + xi{y) — 



(A4) 



y5 + y'' ' " yo + 

This ansatz has been used to obtain the scaling functions shown in Fig. 1 for — 0.5< z <2.0. For |z| outside this interval 
the asymptotic expressions xi{y) and Xs{y) have been used. The constant C2 and the critical exponents /3, 7, S are 
given in Table I, the other parameters needed for this interpolation are collected in Table III. 





ci + da d2 


a b 


yo P 


0(2) 


0.352(30) 0.056 


1.260(3) -1.163(20) 


2.5 3 


0(4) 


0.359(10) -0.025(10) 


1.071(4) -0.866(38) 


5.0 3 



TABLE III: Parameters of the fits to the scaling functions for 0(2) and 0(4). 



APPENDIX B: CHIRAL CONDENSATES FROM CALCULATIONS ON LATTICES WITH TEMPORAL 

EXTENT N-r =4 

In this appendix we present data of our calculations performed with the p4 staggered fermion action on lattices 
with temporal extent Nr = 4 and spatial extent No- = 8, 16 and 32. All calculations have been performed with 2 light 
quarks and a strange quark of mass rhs = 0.065. The improved gauge and fermion actions used for these calculations 
have been described in detail in Ref. [35] . The tables give results of calculations performed at different values of the 
gauge coupling (/3). Results for light (l) and strange (s) quark condensates as well as the disconnected and connected 
contributions to the corresponding susceptibilities are normalized to a single flavor. The last column gives the number 
of trajectories generated for each parameter set. 

In Table VI we give for some selected values of the gauge coupling /? the conversion to a reduced temperature scale. 
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3.3025 


0.0797(5) 


5.63(8) 


6.04(57) 


0.2295(2) 


1.568(2) 


0.89(9) 


16870 


3.3050 


0.0690(4) 


6.47(6) 


5.17(35) 


0.2249(1) 


1.599(1) 


0.69(5) 


21280 


3.3075 


0.0621(7) 


6.94(9) 


4.98(45) 


0.2217(2) 


1.616(2) 


0.63(6) 


6570 


3.3100 


0.0545(7) 


7.45(9) 


4.74(78) 


0.2183(3) 


1.633(2) 


0.72(14) 


7370 






iV^ X Nr 


= 16^ X 4, mia = 0.0016250 








3.2800 


0.1380(7) 


3.77(12) 


2.12(13) 


0.2593(3) 


1.401(2) 


0.44(3) 


21180 


3.2900 


0.1165(7) 


4.47(12) 


3.10(17) 


0.2475(3) 


1.462(2) 


0.55(3) 


27440 


3.3000 


0.0880(10) 


5.74(11) 


5.35(38) 


0.2337(4) 


1.541(3) 


0.76(6) 


40000 


3.3050 


0.0688(10) 


7.11(10) 


5.83(40) 


0.2252(4) 


1.599(3) 


0.82(7) 


42000 


3.3100 


0.0533(11) 


8.14(13) 


4.88(30) 


0.2182(4) 


1.636(3) 


0.63(4) 


24910 


3.3200 


0.0334(7) 


9.05(8) 


2.95(14) 


0.2075(3) 


1.692(2) 


0.51(3) 


25050 


3.3300 


0.0202(4) 


8.16(9) 


1.09(6) 


0.1975(3) 


1.740(2) 


0.39(3) 


14870 






iV^ X Nr 


= 16^ X 4, mia = 0.0032500 








3.2800 


0.1487(4) 


2.84(8) 


1.73(8) 


0.2615(2) 


1.374(3) 


0.44(2) 


40360 


3.2850 


0.1421(9) 


2.93(10) 


1.98(20) 


0.2575(5) 


1.405(5) 


0.49(5) 


13260 


3.2900 


0.1308(5) 


3.20(3) 


2.20(11) 


0.2510(2) 


1.433(2) 


0.50(2) 


45080 


3.2950 


0.1204(9) 


3.56(7) 


2.69(17) 


0.2454(4) 


1.469(4) 


0.60(4) 


19110 


3.3000 


0.1083(6) 


3.88(3) 


3.16(18) 


0.2388(3) 


1.504(2) 


0.64(3) 


41050 


3.3050 


0.0933(7) 


4.43(5) 


3.97(21) 


0.2312(3) 


1.551(3) 


0.76(4) 


39960 


3.3100 


0.0792(7) 


4.82(5) 


4.12(18) 


0.2239(3) 


1.586(3) 


0.77(3) 


42890 


3.3200 


0.0542(6) 


5.88(4) 


3.16(14) 


0.2107(2) 


1.668(2) 


0.66(3) 


44490 


3.3300 


0.0336(3) 


6.25(2) 


1.41(8) 


0.1984(2) 


1.734(1) 


0.46(2) 


39320 







N^ X Nr 


= 32^ X 4, mia = 0.0032500 






3.2800 


0.1488(2) 




1.61(13) 


0.2615(1) 


0.46(3) 


20000 






iV^ X Nr 


= 16^ X 4, mia = 0.0008125 






3.3000 


0.0627(25) 




7.64(94) 


0.2279(9) 


0.71(10) 


6690 






Nj^ X iV, 


= 8=^ X 4, mia = 0.0008125 






3.3000 


0.0452(21) 




4.68(39) 


0.2335(13) 


0.70(5) 


25830 






Af^ X Nr 


= 8=^ X 4, mia = 0.0016250 






3.2800 


0.1141(29) 




5.26(25) 


0.2584(12) 


0.51(4) 


38280 


3.2900 


0.0963(21) 




5.24(17) 


0.2485(9) 


0.61(3) 


40660 


3.3000 


0.0753(35) 




5.32(34) 


0.2380(16) 


0.75(6) 


40100 






Nj^ X iV, 


= 8=^ X 4, mia = 0.0032500 






3.3000 


0.0954(21) 




3.49(13) 


0.2372(10) 


0.63(3) 


30000 



TABLE IV: Light and strange quark condensates {{'4>''p)i.s) for mj/ms < 1/20 and the corresponding disconnected and connected 
parts of the chiral susceptibihties. The last column gives the number of trajectories of half unit length generated for each 
parameter set. 
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/3 




con 

Xi 


dis 
Xl 






xf' 


# traj. 






X Nr 


= 16^ X 4 


mia = 0.0065000 






3.2800 


0.1660(4) 


2.37(1) 


1.11(6) 


0.2661(2) 


1.352(1) 


0.37(2) 


27550 


3.2850 


0.1595(5) 


2.50(2) 


1.14(8) 


0.2619(2) 


1.371(2) 


0.38(2) 


19150 


3.2900 


0.1507(4) 


2.60(1) 


1.37(7) 


0.2562(2) 


1.402(2) 


0.43(2) 


30160 


3.2950 


0.1439(5) 


2.74(2) 


1.77(12) 


0.2519(3) 


1.424(2) 


0.53(3) 


24880 


3.3000 


0.1351(5) 


2.85(1) 


1.94(12) 


0.2464(2) 


1.452(2) 


0.57(3) 


36100 


3.3050 


0.1269(5) 


3.01(4) 


2.23(12) 


0.2414(3) 


1.471(5) 


0.63(3) 


40230 


3.3100 


0.1146(6) 


3.22(4) 


2.74(15) 


0.2343(3) 


1.515(5) 


0.74(4) 


40440 


3.3150 


0.1007(6) 


3.53(3) 


3.09(15) 


0.2262(3) 


1.555(4) 


0.81(4) 


45570 


3.3200 


0.0895(6) 


3.80(2) 


2.53(15) 


0.2197(3) 


1.597(3) 


0.65(4) 


33360 


3.3300 


0.0666(8) 


4.29(2) 


2.10(14) 


0.2061(4) 


1.678(3) 


0.62(4) 


18060 






iV^ X Nr 


= 16''' X 4 


rrna = 0.0130000 






3.2800 


0.1899(3) 


1.98(0) 


0.59(3) 


0.2722(2) 


1.317(1) 


0.28(1) 


20050 


3.2900 


0.1795(4) 


2.07(0) 


0.79(4) 


0.2645(2) 


1.350(1) 


0.34(1) 


21040 


3.3000 


0.1687(5) 


2.18(1) 


1.06(7) 


0.2565(3) 


1.387(3) 


0.45(3) 


18880 


3.3050 


0.1621(4) 


2.27(1) 


1.03(6) 


0.2518(2) 


1.413(4) 


0.43(2) 


32170 


3.3100 


0.1555(6) 


2.31(1) 


1.05(9) 


0.2470(4) 


1.429(2) 


0.44(3) 


16580 


3.3200 


0.1383(6) 


2.55(2) 


1.67(10) 


0.2354(3) 


1.492(4) 


0.64(4) 


28740 


3.3250 


0.1295(4) 


2.72(3) 


1.81(9) 


0.2295(2) 


1.530(13) 


0.68(3) 


54840 


3.3300 


0.1186(5) 


2.81(1) 


1.92(9) 


0.2222(3) 


1.566(2) 


0.72(3) 


50000 






iV^ X Nr 


= 16^ X 4 


rrna = 0.0260000 






3.2800 


0.2256(2) 


1.62(0) 


0.39(2) 


0.2812(1) 


1.267(1) 


0.24(1) 


20310 


3.2900 


0.2180(2) 


1.66(0) 


0.42(2) 


0.2747(2) 


1.290(1) 


0.25(1) 


23950 


3.3000 


0.2105(3) 


1.71(0) 


0.45(2) 


0.2683(2) 


1.314(2) 


0.27(1) 


15330 


3.3050 


0.2054(3) 


1.73(0) 


0.44(3) 


0.2642(2) 


1.331(1) 


0.27(1) 


22550 


3.3100 


0.2010(3) 


1.76(1) 


0.57(4) 


0.2605(2) 


1.346(4) 


0.34(2) 


20170 


3.3200 


0.1916(4) 


1.82(0) 


0.67(4) 


0.2528(3) 


1.378(1) 


0.39(2) 


20030 


3.3300 


0.1807(4) 


1.90(0) 


0.74(6) 


0.2440(3) 


1.417(2) 


0.41(3) 


23380 



TABLE V: Same as Table IV but for the heavier quark masses, mi/ms > 1/10. 





3.2800 3.2900 3.3000 3.3100 3.3200 3.3300 


(T - Tc)/n 


-0.0332 -0.0170 0.0000 0.0181 0.0379 0.0595 



TABLE VI: Relation between (3 and (T - Tc)/Tc for /3c = 3.3000. A shift in /3c of 0.001 corresponds to a shift in (T - rc)/rc 
of about 0.0017. 



